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ABSTRACT 

We present an analysis of the effects of luminosity on the shape of the mid-infrared spectral energy 
distributions (SEDs) of 234 radio-quiet quasars originally presented by Richards et al. In quasars 
without evident dust extinction, the spectrally integrated optical and infrared luminosities are lin- 
early correlated over nearly three decades in luminosity. We find a significant (> 99.99% confidence) 
correlation between the 1.8-8.0/zm spectral index and infrared luminosity that indicates an enhance- 
ment of the mid-infrared continuum with increasing luminosity. Coupled with strong evidence for 
spectral curvature in more luminous quasars, we conclude this trend is likely a manifestation of the 
'near-infrared (3-5/xm) bump' noticed in earlier quasar SED surveys. The strength of this feature 
is indicative of the contribution of emission from the hottest (> 1000 K) dust to the mid-infrared 
spectrum; higher luminosity quasars tend to show more hot dust emission. Finally, the comparable 
distribution of bolometric corrections from the monochromatic 3/im luminosity as well as its lack of 
sensitivity to dust extinction as compared to the standard bolometric correction from ^ sug- 

gests that the former may be a more robust indicator of bolometric quasar luminosity. The close link 
between the power in the mid-infrared and optical and the effect of luminosity on the shape of the 
mid-infrared continuum indicate that considering mid-infrared emission independent of the properties 
of the quasar itself is inadequate for understanding the parsec-scale quasar environment. 
Subject headings: galaxies: active — quasars: general 



1. INTRODUCTION 

Empirically, quasars are remarkable for spectral en- 
ergy distributions (SEDs) with substantial power span- 
ning many decades of frequency from t he radio through 
the hard X-ray (e.g.. lElvis etaflll994h . This SED en- 
codes (though perhaps not unambiguously) the funda- 
mental physical properties of the accreting black hole 
system such as the mass, spin, and accretion rate as well 
as random effects due to orientation; determining the 
mapping from the physical properties of interest to their 
spectral signatures is the ultimate goal of quasar SED 
studies. 

The accretion flow near the black hole generates the 
vast majority of the continuum emission from the near- 
infrared through the hard X-rays, while infrared contin- 
uum light longward of ~ 1 /im is generally attributed 
to a parsec-scale, cold, and dusty region that is heated 
by and reprocesses the higher energy, direct continuum 
and serves to obscure the inner accretion system and 
broad-line region in type 2 (here defined as narrow-line) 
quasars. While the community gcncrically refers to this 
medium as a 'torus', this convenient picture merely ac- 
counts for the empirical evidence for axisymmetric dust 
obscuration without providing a coherent physical de- 
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scription of this structure (see iKonigl fc Kartieiri994l and 
lElitzur fe Shlosmanl 120061 for possible dynamical mod- 
els). Regardless of the specific model, the quasar torus 
might be considered the interface between a quasar and 
its host galaxy, a reservoir for gas to fuel the black hole 
and potentially a site of star formation. Understanding 
the nature of the torus is a primary aim of infrared quasar 
st udies. 

Richards et alj (|2006al hereafter R06) recently pre- 
sented 24/itm through u-band photometric data of the 
largest sample to date of type 1 (broad-line) quasars. The 
resulting 259 SEDs, supplemented with archival VLA, 
GALEX, and ROSAT data as available, span a bolo- 
metric luminosity range of approximately three decades. 
Overall, SEDs of typ e 1 quasars are remarkably simi- 
lar out to z = 6 (e.g.. IShemmer et alj|2006t iHines et"aT] 
12006b I Jiang et alJl200i£ However, in a qualitative com- 
parison of two composite SEDs made from higher and 
lower luminosity objects, R06 noted subtle differences 
indicating that luminosity affected overall SED shapes, 
particularly in the mid-infrared (1.3-8 /im) as shown in 
Figure [TJ Such differences in the mid-infrared spectra 
of quasars with higher mid-infrared fluxes were noted 
previously in the SEDs o f a smaller sample of 3 5 opt i- 
cally selected quasars by lHatziminaoglou et al.1 (2005) , 
and even earlier in an S ED study of 29 Seyfert galaxies 
(|Edelson fc Mall^dfl986h . 

Luminosity is already known to affect some quasar 
emission properties. For example, the Baldwin Effect 
characterizes the decreasing equivalent widths of ul- 
traviolet emission lines with increasing ultraviolet con- 
tinuum luminosity (e.g ., iBaldwml Il977t Richard s et al.l 
2002; Sh ang et al]|2003h . In addition, the maximum ve- 
locity of outflows seen in C IV absorption increases with 
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ultraviolet luminosity (Laor fc Brandt! l2002fh Across 
broader spectral regions, the X-ray-to-optical flux ra- 
tio is known t o decrease with increasing ultraviolet lu- 
mino sity fe.g.. lAvni fc Tananbaum|[l986t IStrateva et al.1 
2005). Finally, and perhaps most relevant for the cur- 
rent study, the ratio of type 1 to type 2 quasars is 
found to increase with luminosity, as fou nd in X-ray 
fe.g., lUeda et all 12003); lAkvlas et al.ll2006D and narrow 
emission-line surveys ( Hao et al.ll2005af) . Overall, these 
empirical relationships indicate that luminosity affects 
the quasar system in some fashion. 

In this paper, we investigate potential SED differences 
of type 1 quasars as a function of luminosity, focusing in 
particular on the rest-frame mid-infrared and optical/UV 
emission with the goal of understanding how the quasar 
affects its immediate environment. Throughout we as- 
sume A-CDM cosmolog y with £7m = 0-3, = .7, and 
Ho = TOkms-iMpc" 1 ilSpergel et al.ll2003ll2006h . 

2. ANALYSIS AND RESULTS 

We take as the starting point for our analysis the pho- 
tometric SED data from Tables 1 and 2 of R06, fo- 
cusing on the spectral regions with complete coverage: 
the observed-frame optical and mid-infrared. The opti- 
cal coverage is fro m the Sloan Di gital Sky Survey ugriz 
photometric data (|York et al.ll2000f) corrected for Galac- 
tic extinction. The mid-infrared data are from Spitzer 
IRAC 3.6, 4.5, 5.8, and 8.0/^m (full widths at half maxima 
[FWHM] of the spectral responses are 0.7, 1.0, 1.4, and 
2.9/Ltm, respectively) and MIPS 24/^m (FWHM=5.3 / um) 
imaging. There has been no attempt to remove host 
galaxy light, which is likely to contribute most signifi- 
cantly to the SED at rest-frame ~ 1/im. 

R06 included all SDSS quasars with publicly avail- 
able Spitz er coverage as o f the D ata Release 3 quasar 
catalog of ISchneider et all (|2005h . The SDSS is quite 
complete to all but t he reddest quasars to i = 19.1 
(iRichards et alJl2006blk all SDSS quasars were detected 
by Spitzer. From the entire R06 sample of 259 quasars, 
those that are known to be radio luminous (log(Lp{ ) > 32 
ergs -1 Hz -1 at rest-frame 20 cm assuming a„,R = —1) or 
have no 24 ^m photometry have been eliminated. Note 
that the majority of R06 quasars (229/259) have only ra- 
dio upper limits. Of these 229, 99 do not have sensitive 
enough flux limits to rule out high radio luminosity. The 
trimmed sample comprises 234 quasars with z — 0.14- 
5.21, z = 1.54 ±0.92, z = 16.8-20.4, i = 18.9 ± 0.7, 
Mi =-22.0 to -28.2, and M, = -25.6 ± 1.6. 

Throughout this section, when we refer to a linear fit, 
we are using the ordinary least-squares bisector method 
to determine the slope and inter cept of a line fitte d to 
(usually logarithmic) data values. Ilsobe et al.l (|1990f ) rec- 
ommend this method for determining the functional rela- 
tionship between ordinate and abscissa values when the 
errors are not well defined or insignificant. Errors quoted 
for the slopes are the standard deviations of the best- 
fitting values. 

2.1. Relative Integrated Infrared and Optical 
Luminosities 

As a first step, we investigated whether the relative 
power in the infrared vs. the optical varies with lumi- 
nosity. As mid-infrared emission and optical photons are 
generated in physically distinct regions, there is no rea- 



son a priori why the relative power should (or should 
not) be affected by luminosity. Figure [2^ shows the log- 
arithm of the ratio of infrared and optical luminosities vs. 
the logarithm of the optical luminosity; the quantity /i r - 
l opt may be considered a parameterization of the repro- 
cessed to direct accretion-powered luminosity. An alter- 
nate depiction of the vs. l op t is presented 
in Figure [2p. We use k T and l opt to indicate log(L op t) 
and log(ii r ), respectively, with units of erg s _1 ; Both 
values are taken from Table 1 of R06. The first, k Y , 
is the luminosity integrated from log(i/) = 12.48-14.48 
Hz (l-100^im); l opt is integrated from log(z^) = 14.48- 
15.48 Hz (0.1-l^m). Power-law interpolation was used 
on the monochromatic luminosities calculated from the 
tabulated flux densities. 

From visual inspection, there is a clear trend with in- 
creasing optical luminosity for quasars to have relatively 
more power in the optical, and a linear fit gives a slope 
of —0.51 ± 0.06 (la uncertainty). Non-parametric sta- 
tistical tests of the bivariate data measuring Spearman's 
p=-0.288 and Kendall's r=-0.194 indicate that the prob- 
ability of obtaining these values if the parameters are not 
correlated is < 1 x 10 -5 ; the statistical values and prob- 
abilities are listed in Table [1] 

Given that l op t characterizes direct emission from the 
accretion disk, it may be considered in some sense a 
truer measure of overall quasar power than the indi- 
rect and reprocessed infrared light. A straightforward 
interpretation of this significant trend would be that a 
smaller fraction of the accretion-powered luminosity is 
reprocessed and emitted in the mid-infrared in more lu- 
minous quasars. However, optical photons are subject 
to dust reddening and extinction. The integrated optical 
luminosity measured by the observer would be reduced 
by extinction while h r -l pt would increase in the same 
sense as the correlation. Therefore, if dust extinction 
is prevalent, it could induce the correlation seen in Fig- 
ure [2^. As the small grains responsible for moderate 
optical extinction will not affect the infrared power (and 
may not in any case be along the same line of sight), l- 1T 
may be a more robust measure of the total quasar power. 
Indeed, we perform the same statistical tests on /j r -/ pt 
vs. l lr and find that these two quantities are not corre- 
lated. The scatter about the mean is also roughly the 
same across almost three decades in Z; r . 

2.2. Comparison of Optically Blue and Red Quasars 

The values of Zi r -Z pt evident in Figure^ span roughly 
one decade from -0.5 to 0.5. As the apparent trend of 
'ir-^opt vs. l ov t plotted in Figured may be induced by 
extinction, we investigate if extinction could also con- 
tribute to the width of the spread in k r -l pt- To simply 
parameterize the color of the optical continuum, a power- 
law model (In oc v a "< as *') was fit to the SDSS photometric 
data between rest-frame 1200-5000 A. This region of the 
spectrum is typically dominated by the quasar, and by 
inspection of the continuum fits compared to the data 
(2-5 points in this region); a power-law model is appro- 
priate for this purpose. For illustration, three sample 
SEDs are shown in Figure [3l From analysis of the SDSS 
composite quasar spectrum, the host galaxy contribu- 
tion reddens the optical c ontinuum slope for A > 5000 A 
(jVanden Berk et al.|[200lh . Wavelengths < 1200 A are to 
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be avoided as the intervening Lya forest reduces the flux 
significantly. 

A histogram of a^opt values is shown in Figure^. At 
the blue (more positive) end, the distribution is roughly 
gaussian; a long tail extends to the red (negative) end. 
The qualitative shape of this a„ i0 pt distribution is con- 
sistent with that of A(g — i), the SDSS g — i color 
minus the mo de of the g — i colo r of quasars at the 
same redshift (jRichards et al.l [2003h . The very red tail 
in the A(g — i) distribution is associated with contin- 
uum dust reddening, though intrinsically red quasars 
may contribute somewhat. To illustrate the effect of op- 
tical extinction on Figures [2^, a reddening vector has 
been drawn. The vector reflects the decrease in l opt 
of A g = 0.1 magnitudes of extinction, assuming Small 
Magellanic C loud (SMC) reddenin g (as appropriate for 
quasars; e.g.. iHopkins et aLll2004h of an a„ opt = —0.49 
optical continuum integrated from 1000-10,000 A (fol- 
lowing R06). 

We use the median a VjOp t=—0A9 to divide the sample 
into 'optically red' (a„, pt < —0.49) and 'optically blue' 
(ovopt > —0.49) groups. A comparison of the histograms 
of iir-^opt for the groups (Figure [4p) shows an offset. A 
student's t-test of the two populations indicates that the 
probability that they have the same mean is < 1 x 10~ 6 . 
This mean reduction in l- 1T — l opt for the red quasars is 
consistent with the expectation from extinction. In fact, 
within the optically red quasars, the subset of clearly 
dust-reddened {ot v , -pt< —0.70) quasars show the most 
extreme values of l a -l op t- 

Therefore, we conclude that optical extinction makes 
a significant contribution to the width in the li t -l op t dis- 
tribution, and likely induces the correlation of Zi r — Z opt vs. 
Z opt . In fact, this correlation disappears when only blue 
quasars are considered (see Table [lj . For optically blue 
quasars, arguably the least likely to suffer from dust red- 
dening, the data are consistent with a constant ratio of 
infrared to optical luminosity over approximately three 
decades in luminosity: L 1T oc L opt (see Figure EfcO- 

2.3. Dependence of the Mid-Infrared SED on 
Luminosity 

A close-up of the composite SEDs from R06 in the spec- 
tral region from 1-12 /im reveals distinctions between the 
luminous and faint quasar SEDs (Figure [T]) . When the 
SEDs are all normalized to 1.3 /im, the luminous SED 
shows enhanced flux in the mid-infrared spectral region 
relative to the mean and faint composites. The three 
SEDs then almost meet again for A >10/im. Note that 
these composites were constructed from broad-band pho- 
tometry, and therefore do not reflect interesting mid- 
infrared features such as strong 10/im silicat e emission 
recentl y seen in several qu asar spectra (e.g., lHao et aTl 
I2005bt ISiebenmorgen et~aT]|2005l ). 

To investigate if these differences arc evident trends in 
all of the data rather than just in the composites, power- 
law fits have been performed to the Spitzer data in the 
rest-frame 1.8-8 /im region. As available, 15/im ISO and 
JHK photometry were also incorporated. This spectral 
region was chosen to maximize the photometric cover- 
age while reducing to the extent possible contamination 
from the host galaxy (which peaks in the near-infrared) 
and silicate emission (at ~ 10/im). Even in logarith- 
mic units, curvature is evident in this spectral region 



(Figure [TJ, and so the spectral index a v ^ (l v oc v aiJ <") 
should only be considered a simple parameterization of 
the data. Sample SEDs and the corresponding fits are 
shown in Figure [3l The value of a v ^ t for eight quasars 
with z > 3.5 could not be measured because the MIPS 
24 /im data point was the only one in the bandpass. 

A histogram of the distribution of a u ^ values is shown 
in Figure[4^. For reference, the mean and standard devi- 
ation of the OLujx values, —1.20 ±0.34, are con sistent with 
the 1 -10/im spectral indices measured by lHaas et aTl 
(2003) of —1.3 ± 0.3 for 64 Palomar-Green quasars ob- 
served with ISO. 

As anticipated from the composite spectral results, the 
distribution of a u ^ v vs. l lr shows a clear trend with in- 
creasing l[ T whereby the more luminous quasars tend to 
have steeper (more negative) spectral indices in the mid- 
infrared; this is evident in Figure [5j This apparent cor- 
relation is supported with the bivariate Spearman and 
Kendall tests, which give probabilities that a v ^ is not 
correlated with l ir of 6.9 x 10" 5 and 2.7 x 10~ 5 , respec- 
tively. (We use kr instead of l opt to avoid the problems 
of underestimating the luminosity due to extinction.) 

Next, we investigate if quasars are also brighter at 
1.8-8 /xm relative to the optical continuum region. We 
choose to measure the monochromatic luminosities at 
3 um and 5000 A (U„ m and / ° in logarithmic units) 

V o/im 5000A ° ' 

to represent the mid-infrared and optical, respectively. 
Three /im is near the middle of the 1.8-8.0 /im bandpass, 
has good photometric coverage, and should have mini- 
mal contribution from even a star-forming host galaxy; 
5000 A marks the red end of the quasar-dominated op- 
tical/UV continuum region. As seen from the results in 
Table [1] l^um-l "does not change significantly as h ura 

5000A ^ 

increases. If only optically blue quasars are considered 
(to reduce the potential effects of extinction), there is 
still no significant evidence for any change in lR,, m -l 

to _ . 5000A 

for increasing Iz^m- The steepening of the mid- infrared 
continuum evident from vs. k r is much clearer evi- 
dence of mid-infrared changes. 

In the course of the previous tests, we visually ex- 
amined each of the best-fitting optical and mid-infrared 
power-law fits in comparison with all of the SDSS and 
Spitzer photometry (cf., Figure [3]) ■ Perhaps surprisingly, 
up to z r-j 1.2, extrapolating a simple power-law fit from 
the IRAC photometry predicts quite well L v in the MIPS 
24/tm bandpass. From z = 1.2-2, however, the extrapo- 
lation typically overpredicts L u in the MIPS 24/iin band- 
pass. (Beyond 2 = 2, the MIPS datapoint is included in 
the power-law fit.) This systematic effect is not because 
of a contribution of significant 10/im flux from silicate 
emission, which would cause the opposite effect. Instead, 
the deficit of 24/im emission is evidence for curvature in 
the mid-infrared spectrum. This is consistent with the 
luminous and faint composite SEDs, which when normal- 
ized at 1.3/xm almost meet again near 10/im (see Fig. Q] 
and R06). 

To characterize spectral curvature, we calculate the 
monochromatic 8/im luminosity, -Ls/xm, in two ways. For 
the first, Ls t a, the best-fitting mid-infared power-law 
spectral index and normalization are used to extrapo- 
late to rest-frame 8/im. For the second, £s,SEDi photo- 
metric datapoints within 0.3 dex of rest-frame 8/im are 
used to normalize the R06 composite SED to measure 
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the 8/im luminosity. The lowest panel in Figure [3] illus- 
trates the discrepancy between Ls,sed and Ls, a seen in 
some objects. Quasars without at least three data points 
in the 1.8-8,0/im regime and photometry within 0.3 dex 
of 8/im were eliminated from this analysis to mitigate 
redshift effects, the resulting sample has 1 46 objects. As 
seen in Figure[5J the strong decrease in L s ^sed/Ls,c with 
increasing L lr is indicative of significant curvature in the 
spectrum (see Tabled]). 

Overall, we find strong evidence for significant differ- 
ences in the mid-infrared spectral regime as a function 
of luminosity. Specifically, more luminous quasars show 
steeper (more negative) values of a„,; r . Extrapolating 
a. v ; xt to rest-frame 8/im also overpredicts the 8/im lumi- 
nosity in luminous quasars, thus demonstrating spectral 
curvature. Both these results are consistent with the 
expectation from the luminous composite shown in Fig- 
ure [TJ 

2.4. Investigating Potential Systematic Redshift Effects 

Quasar spectral studies have historically been plagued 
with difficulties in distinguishing betwe en luminosity and 
redshift effects (see lSteffen et al.ll2006l for a recent exam- 
ple) particularly as spectral redshifting causes distinct 
spectral regions to be sampled in a given observed-frame 
bandpass. 

In this study, the evident mid-infrared spectral cur- 
vature and the gap in the Spitzer photometric coverage 
between observed-frame 8 and 24/im could induce appar- 
ent systematic differences in the mid-infrared spectra as 
a result of redshift. To examine this potential problem, 
we focused on a narrow redshift slice in our sample, and 
reran the statistical tests of a u ; lr vs. l a . We choose the 
redshift range from z — 1.0-1.5 because there is still a 
significant (> 1 dex) infrared luminosity range and at 
least 3 IRAC bands sample the rest-frame wavelength 
range of interest. This redshift subsamplc comprises 55 
quasars. Even for this significantly smaller sample, the 
correlation of a Vi i T with l\ T remains statistically signifi- 
cant, with both Spearman's and Kendall's tests giving 
< 1 x 10~ 4 probabilities of no correlation. Furthermore, 
if the effect were entirely induced by bandpass effects, 
ov,ir would be expected to correlate significantly with 
redshift, which is not the case (see Table [TJ. More even 
photometric sampling is likely needed to better constrain 
the shape of the correlation and certainly to investigate 
spectral curvature in greater detail. 

2.5. Comparing Optical and Infrared Spectral Indices 

In the standard understanding of the infrared emission 
in quasars, cold, dusty material on parsec scales from the 
supcrmassive black hole is illuminated by and reprocesses 
accretion disk radiation. The torus geometry is typically 
described as some form of flattened disk-like configura- 
tion, encompassing toroids, open cones, and flared disks. 
At mid-infrared (1-10/iin) wavelengths, the continuum 
spectrum is dominated by thermal dust emission, and 
as long as the highest dust temperature is fixed (by dust 
sublimation, for example), the SED is not sensitive to the 
spectral shape of the illumi nating continuum reg ardless 
of the assumed geometry (jlvezic fc Elitzurlll997 ). Only 
at near-infrared (< 3/im) wavelengths is the illuminat- 
ing spectral shape expected to affect the SED as a re- 
sult of scattering of direct continuum photons (Nenkova 



et al., in prep.). As seen in Figure [3 from these data 
we find no evidence for any correlation between a„ !0 pt 
and a^ir, confirming theoretical expectation. Even when 
considering just the optically blue quasars to reduce the 
potential masking effects of optical extinction, the two 
parameters show no relation. Therefore, there is no evi- 
dence that the mid-infrared spectral changes are induced 
by the shape of the illuminating continuum; this is also 
consistent with the lack of notable differences between 
the composite SEDs constructed by R06 from optically 
red and optically blue quasars. 

2.6. Evaluating Mid-Infrared and Optical Bolometric 
Corrections 

Typically, quasars only have photometric data in a few 
bandpasses while the quantity of interest (e.g., for esti- 
mating black hole masses and Eddington accretion rates) 
is the integrated bolometric luminosity. Therefore, find- 
ing a single point in the SED from which a robust es- 
timate of the bolometric luminosity can be made would 
be useful. Given the evidence discussed in $'2.1\ for the 
potentially significant effects of dust extinction on l op t, 
we investigate whether a mid-infrared point would be 
more robust. Using the 3/im monochromatic luminosity, 
vLsfan, we calculated both the overall (to Lboi) bolo- 
metric corrections (BCs) and the infrared (to L- u ) BCs. 
The value of Lboi is integrated from 100/xm to 10 keV; 
where photometry do not exist to measure this quan- 
tity, the composite SED was used to fill out the wave- 
length coverage. These are also compared to the typical, 
5100A BCs as compiled in R06. Histo grams of the three 
BC distributions are shown in Figure [51 The 3/Ltm BC, 
8.59±3.30 is notably quite comparable to the 5100A BC, 
10.47±4.14, with fractional standard deviations of 38% 
and 40%, respectively. The extremely tight and symmet- 
ric 3/im infrared BC indicates that vL^^ m is quite a good 
proxy for L 1T . 

3. DISCUSSION 

The distribution of l lr vs. l op t as plotted in Figure[2}3 is 
notably quite tight. As expected in a paradigm where the 
direct optical/UV continuum emission powers thermal 
emission in the infrared, optical and infrared luminos- 
ity increase in tandem. While the overall sample of 234 
quasars gives L ir cx x°- 94±0 02 , considering only the blue 
quasars (the least likely to be dust-extincted) gives a lin- 
ear relationship. This constant ratio of optical to infrared 
luminosities over ~ 3 decades in luminosity suggests that 
the dust mass (estimated from the infrared luminosity) 
radiating in the regime probed by Spitzer also increases 
approximately linearly with optical luminosity (assum- 
ing similar grain properties and distributions). At the 
same time, the inner radius of the dust-emitting region 
increases as V-^uv as the dust sublimation radius moves 
out. Thus the inner wall of the torus at any luminosity is 
expected to have a similar temperature. It is challenging 
to extrapolate this point to inferences about the total 
infared emitting volume, as for optically thick emission 
(expected for A < 10/im), the observed luminosity does 
not depend only on volume, but also on, for example, 
dust covering fraction, self-shadowing, and viewing an- 
gle. Therefore, it is not clear how our results are related 
to the interpretation that torus geometry must change 
with quasar power, as understood from the increasing 
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type 1 /type 2 ratio with luminosity found e mpirically in 
X-ray (|Ueda et al]l2003l: lAkvlas et al.ll2006l ) and narrow 
emission-line surveys (jHao et al J l2005alh 

Indications that quasar luminosity is affecting the torus 
structure are the significant spectral changes in the mid- 
infrared. Specifically, we seen an enhancement of mid- 
infrared emission in more luminous quasars manifested as 
a steepening of the 1.8-8.0/im continuum. As expected 
from the composite SEDs, this enhancement is better de- 
scribed as a bump, as an extrapolation from a„.i r over- 
predicts the luminosity at ~ 8/im in luminous quasars. In 
current static torus models, the shape of the mid-infrared 
continuum shortward of the 10/im silicate feature is sen- 
sitive to numerous physical inputs , inclu ding (1) orien- 
tation effects (e.g.. iNenkova et ail I2002D. (2) the dust 
grain properties (e.g., see Fig. 6 of iSchartmann et alj 
120051) . (3) the dumpiness and num ber of emitting clouds 
( Dullemond fc van Bemmel 120051 : Nenkova et al., in 
prep.), an d (4) the openin g angle of the torus (e.g., see 
Fig. 20 of lFritz et a.l.ll200 6h . Our focus on type 1 quasars 
has reduced the effect of (1) to the extent possible. Sim- 
ilarly, given that the inner wall of the torus is set by the 
graphite sublimation temperature to be T ~ 1500 K, 
a luminosity dependence of dust grain size, composi- 
tion, and distribution is not expected, though detailed 
infrared spectroscopy is required to examine this issue 
empirically. Regarding (3), current clumpy torus models 
do not generate explicit predictions for any luminosity 
dependence of the mid-infrared continuum. Therefore, 
option (4) seems plausible given the other evidence for 
a luminosity-dependent covering fraction of the torus. 
However, we are not aware of explicit predictions in the 
literature for how a changing opening angle would affect 
the mid-infrared continuum that are consistent with our 
results. 

Another possible explanation for the enhanced mid- 
infrared spectra of luminous quasars is that a hotter dust 
component may be contributing more significantly to 
their SEDs than in lower luminosity objects. A striking 
exa mple of this type of component w as rec ently presented 
by iRodrfguez-Ardila fc Mazzalavl ((2006) in the near- 
infrared spectrum of the narrow-line Seyfert 1 galaxy, 
Mrk 1239; they interpreted it as thermal emission from 
T ~ 1200 K graphite grains. This type of spectral fea- 
ture, the so-cal led 'near-infrared bum p', has been noted 
for years fe.g.. iHvland fc Allen! [1982T ) . and it s physical 
origin was interpreted in the same fashion by iBarvainlsl 
(1987). Similar structure is also evident in the mean 
radio-quiet composite spectrum of Elvis et al. (1994; 
though not, interestingly, in the radio- loud composite 
- see t heir Fig. 10), while th e more recent combined 
SEDs of lHatzimin aoglou et al.l (|2005h also show some 
hint of the near-infared b ump (see their Fig. 9). In fact, 
lEdelson fe Malkanl (|1986h (who modeled the feature as 
a parabola peaked at 5.2/im in log(/„) vs. log(i/) units) 
found its strength to be correlated with luminosity in 
their study of 29 quasar SEDs. Given the IRAC wave- 
length sampling, at higher redshifts, the data are probing 
the shorter wavelengths of the mid-infrared bandpass, 
and the presence of a 'hot bump' peaking between 3- 
5/zm would manifest itself as a steepening in a^h - The 
overprediction of Ls/im from extrapolations from the mid- 
infrared power law further supports this interpretation. 
This conclusion is not necessarily inconsistent with a 're- 



ceding torus' paradigm whereby m ore luminous qu asars 
have tori covering less of the sky (|Lawrencdll991[ ) as a 
flatter dusty structure would be less likely to block lines 
of sight to the hottest dust. 

Finally, we also mention that polyaromatic hydrocar- 
bon (PAH) emission is a possible contributor to the mid- 
infrared enhancement in more luminous quasars. We 
consider this unlikely however, as the strongest PAH fea- 
tures (at 6.2/mi and 7.7/mi) are not sampled by the IRAC 
data beyond z ~ 0.3 and do not enter the MIPS 24/im 
bandpass until z ~ 3. Other PAH features in the 2-5/xm 
range are much weaker in star-form ing regions and lab- 
oratory spectra (|Peeters et al.l 12004). though they have 
not been well-studied in luminous quasar spectra. 

The gross similarity of quasar mid-infrare d through op- 
tical S EDs is consistent with the finding of iBrown et al.l 
(2006) that both mid-infrared and optical type 1 quasar 
selection give similar space densities and luminosity func- 
tions for z = 1-5. Therefore, though optical dust extinc- 
tion may induce a false correlation as in Figured, it is 
unlikely to distort strongly the study of type 1 quasars. 
Similar to our findings, a lack of variation of the ratio of 
far-infrared to optical luminosity as a function o f bolo - 
metric luminosity was found by lAndreani et al.l (|1999h 
for a sample of 120 quasars (with many far-infrared up- 
per limits). However, we caution on the use of far- 
infrared (> 20/im) luminosities to characterize quasars 
as star formation in the host galaxy can contribute sig- 
nificantly to or even dominate the power in this regime. 
For example, from the 1-100/im spectral modeling of 36 
type 1 quasars with luminosities matched t o our sam- 
ple (L quasar > 3 x 10 44 erg s _1 ) presented in iFritz et al.l 
(2006), the range of modeled quasar contributions to 
the total 5-1000/im luminosities ranged from 9-100% 
with a mean of 62 ± 26%. The correlation between the 
far -infrared luminos i ties an d 7.7/im PAH emission seen 
by ISchweitzer et al.l (|2006l ) in a sample of quasars and 
ultra-luminous infrared galaxies further supports star- 
formation as a significant (and variable) source of far- 
infrared power in luminous active galaxies. However, in 
typical torus+star formation model fits, the quasar domi- 
nates from ~ 2 to 10/im fe.g-. lvan Bemmel fc Dullemondl 
2003), even considering some contribution from PAH 
emission. As the integrated L lr values used in this work 
were derived from IRAC+24/ini MIPS data, they may 
more accurately represent the quasar power than infrared 
luminosities estimated from longer wavelength data. 

4. SUMMARY AND CONCLUSIONS 

We have analyzed 234 radio-quiet quasar SEDs from 
the sample of R06 with the goal of investigating possible 
mid-infrared SED changes as a function of luminosity. 
While we find quasar SEDs overall to be remarkably sta- 
ble with luminosity in this wavelength regime, there are 
some notable differences. We detail the following results: 

1. Over three decades in infrared luminosity, the mean 
ratio and scatter of Li r /L opt vs. L lr in quasar SEDs 
is constant. Some part of the spread in the dis- 
tribution of L; r /L pt is incurred from optical dust 
extinction. Given that dust extinction can signif- 
icantly reduce i pt, the mid-infrared luminosity 
may be a more robust indicator of bolometric lu- 
minosity. 
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2. The steepening of the 1.8-8.0/zm spectral index, 
cki/,i r , with luminosity coupled with the spectral 
curvature evident from L^^m measurements in- 
dicate that more luminous quasars show a mid- 
infrared enhancement consistent with a 3-5/im 
bump. This result confirms the qualitative differ- 
ences noted by R06 between luminous and faint 
composite quasar SEDs. As noted in previous stud- 
ies, this feature can be attributed to an increasing 
contribution from hot (T > 1000 K) d ust with in- 
creasing quasar power (|Barvainislll987l ). 

3. The overall bolometric correction from the 
monochromatic 3/xm luminosity, 8.59±3.30, has 
a comparable dispersion to the typical bolomet- 
ric correction from vL ° without suffering from 

5100A _ ° 

dust extinction. The infrared bolometric correc- 
tion from vL-s^ to L„, 3.44 ± 1.68, shows a very 
tight and symmetric distribution, and is therefore 
a good measure of mid-infrared power. 



Static models for the torus are ultimately unsatisfy- 
ing; our results support the picture where quasars are 
shaping the environment that emits in the mid-infrared 
- the 'torus' is not simply a passive reprocessor of 
UV/optical emission. The dynamical wind paradigm, 
where magnetohydrodynamic and radiation pressure lift 
and accelerate a dusty wind aw ay from the central en- 
gine is therefore more appealing (IKonigl fe Kartidll994l : 
lEverett et all 120021 : lElitzur fc Shlosmanl l2006h . and ra- 
diation pressure naturally has a luminosity dependence. 
Detailed predictions of mid-infrared spectra within this 
framework are a necessary next step, and future studies 
of 1-10/im infrared spectra of type 1 quasars matched 
in redshift with a range of luminosities will illuminate 
the nature of the spectral differences we find. Compar- 
isons with theoretical models of such data will enable 
constraints on physical properties of interest such as the 
temperature, structure, and radius of the dusty region 
to aid in understanding the role of luminosity in shaping 
the quasar environment. 



Though the first two results point strongly to a change 
in the geometry of the mid-infrared emitting region with 
luminosity, we know of no predictions from current torus 
models that are consistent with these findings. The no- 
table 1.8-8.0/im spectral changes we have found therefore 
furnish strong constraints for future theoretical modeling 
of mid-infrared torus emission. 
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TABLE 1 

Results from Non-Parametric Bivariate Statistical Tests 



Variables Spearman Kendall 



Independent/Dependent (Sample) a 


P 


Prob. b 


T 


Prob. b 


«opt/('ir-«opt) (234) 


-0.288 


7.52e-6 


— 0.194 


9.66e-6 


lopt/ihr-lopt) (blue; 117) 


-0.077 


4.09e-l 


— 0.052 


4.08e-l 


lo P t/(hr-lopt) (red; 117) 


-0.382 


2.18c-5 


-0.260 


3.24e-5 


W('ir-Zopt) (234) 


0.039 


5.51e-l 


0.032 


4.66e-l 


(kr-lopt)/ttv,opt (234) 


-0.422 


< le-7 


-0.292 


< le-7 


l iT /a Vtil (226) 


-0.262 


6.89e-5 


-0.188 


2.67e-5 


h T /a v ir (z = 1.0-1.5; 55) 


-0.508 


7.48c-5 


-0.375 


5.36e-5 


z/a„,i'r (226) 


-0.058 


3.88e-l 


-0.061 


1.74e-l 


av,opt/av,ir (226) 


0.020 


7.64e-l 


0.013 


7.67c-l 


av, op t/av,ir (blue; 113) 
^m/(V m -« 5000A (234) 


0.126 


1.84e-l 


0.088 


1.68e-l 


0.074 


2.62e-l 


0.047 


2.85e-l 


hum/ikum-l i> (blue; 117) 

«ir/(«8,SED-i 8 ,a) (146) 


0.235 


1.07c-2 


0.159 


l.lle-2 


-0.337 


3.17c-05 


-0.224 


5.90c-05 



a The sample labels and numbers of data points are given in parentheses. 'Blue' and 'red' refer to a„ l0p t > —0.49 and ct„ l0 pt < —0.49, 
respectively. 13 The two-sided probability that the given variables are not correlated. 
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Fig. 1. — Plots of three SED composites from R06 made from the entire sample (solid curve), the luminous subset (log(L op t) > 46.02 
ergs -1 ; dashed curve), and the faint subset (log(L opt ) < 46.02 ergs -1 ; dot-dashed curve) in units of (left) vL v , to emphasize the 
difference in spectral shapes, and (right) L v , to show the spectral region as we model it. The composite SEDs have been normalized 
to match at 1.3/im. The SED resolution is a convolution of the IRAC (FWHM=0.7, 1.0, 1.4, and 2.9^tm at 3.6, 4.5, 5.8, and 8.0/xm, 
respectively) and 24/xm MIPS (FWHM=5.3/im) filter functions with the redshift distribution. Statistical subtraction of the host galaxy 
light is included; for detailed information on the composite construction see R06. For reference, 8 and 1.8/im are indicated with vertical 
dotted lines. 




Fig. 2. — (Left) Plot of log(Li r /L op t) vs. log(L op t) for the 234 SED quasars. The solid line represents a linear fit to the plotted data; 
the corresponding slope and standard deviation are given in the figure. The labelled vector indicates the magnitude and direction of the 
effect of 0.1 magnitudes of g-band extinction for SMC reddening of a power-law continuum with the median a„ l0p t= —0.49. (Right) Plot 
of log(Li r ) vs. log(L op t) for the same sample. The optically bluest and reddest thirds are indicated with blue crosses and red diamonds, 
respectively. Linear fits to all the data as well as the bluest and reddest thirds are drawn in solid lines; the best-fitting slopes are marked 
on the figure. While the lefthand panel indicates a significant correlation between the ratio of optical to infrared luminosity and the optical 
luminosity, an examination of the optically blue data imply that the data are consistent with L; r oc L opt if extinction is taken into account. 
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Fig. 3. — Sample SEDs illustrating the data analysis process; the abscissa is rest-frame frequency. Da ta are labeled with crosses, open 
thick diamonds indicate the computed quantities (from left to right) of Lg, Q , Lg^, and L ° (see i|2,3l for definitions). The thin dashed 
curve shows the mean composite SED from R06 normalized to the photometric data within 0.3 dex of rest-frame 8/xm; the thin open 
diamond represents -Lg^ED (overlapping with Lg jC( in the top two panels). Red dashed and blue dotted vertical lines mark the boundaries 
of the mid-infrared and UV/optical fitting regimes, respectively. Solid lines are the best-fitting power-law models to the mid-infrared and 
optical data; they have been extrapolated to cover all of the data. The top panel is from a quasar with data closely matching the composite 
SED. The middle example shows a red (likely from dust reddening - note the decreased L op t) optical continuum, while the bottom panel 
indicates a measurable discrepancy between Lg^ED arltl ^8, a indicating mid-infrared spectral curvature. 
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Fig. 4. — (Left:) Histograms of ctv.opt (solid) and ct, y j r (dot-dashed) for the SED sample. The ol u values are determined from power-law 
fits {I u cc v av ) to the photometric data; bluer continua are to the right. The median value for a^opt of —0.49 is indicated with a solid 
vertical line. For reference, av, QD t=— 0. 44, the best-fit value to the continuum of the composite SDSS quasar spectrum is also plotted (dotted 
vertical line; Vandcn Berk ct al. 2001). (Right:) Histograms of the distribution of log(L; r /Z pt) for the optically blue (av,opt> —0.49; 
solid), red (av j0p t< —0.49; dashed) quasar populations. Quasars that are very likely to be dust-reddened (a l /, pt< —0.70; dotted) are also 
shown. The difference between these histograms indicates that at least some of the spread in the overall log(Li r /L op t) distribution is an 
artifact of optical extinction that reduces L Q pt ■ 
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Fig. 5. — The infrared spectral index, a Vi i T , vs. log(Li r ) for the SED sample (solid dots); the subset of the SED sample with z = 1.0—1.5 
is also marked with crosses. The parameter a v ^ is determined from power-law fits (l v oc v a "< ir ) to the Spitzer photometry between 
rest-frame 1.8— 8.0 /im; the eight quasars at z > 3.5 (plotted at a v ; r =— 3.0) only had one photometric data point in the bandpass; they 
are not included in the analysis. For reference, a horizontal dashed line marks the median a VtlI value; for higher values of L; r a larger 
fraction of quasars lie below the median. The median uncertainty in a v j r from fitting a power-law model to the photometry is indicated 
as a vertical errorbar on a large, filled circle in the upper right corner of the plot. 



Luminosity Effects on Mid-IR Quasar SEDs 



11 



CO 



O 
Ld 

CO 

cd 



CP 

o 



0.4 



0.2 - 



-0.0 ; 



— ' -0.2 - 



-0.4 - 



-0.6 



• • ? * • • * • • 



—I I I I I I I I I I I— 



—I I I I I I I I l_ 



44.0 44.5 45.0 45.5 46.0 46.5 47.0 47.5 
log(U r ) (erg s" 1 ) 



Fig. 6. — The logarithmic ratio of 1/8, sed determined from the R06 mean composite SED normalized to photometry within 0.3 

dex of rest-frame 8/im) to Lg^ a (Ls/jm calculated by extrapolating to 8/^m from the mid-infrared power-law model) versus log(Li r ). For 
reference, a horizontal dashed lines marks a ratio of unity. The strong decrease in -Ls.SEd/ L»,a with increasing L; r indicates increasing 
spectral curvature at higher luminosities. 
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Fig. 7. — The infrared spectral index, av,ir, versus the optical spectral index, a^opt, for the entire SED sample. Quasars with only one 
photometric data point in the 1.8— 8.0/xm bandpass are plotted at «„,ir = —3.0. The symbol sizes have been scaled to indicate the relative 
Li r of each quasar; the larger symbols indicate more luminous quasars. For clarity, quasars with log(Li r )< 46.0erg s _1 (the median of the 
L; r distribution) have been colored green; more luminous quasars are blue. The cross in the lower left corner of the plot indicates the size 
of the median errors for the fitted values of a v ; r and ct^opt- 
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Fig. 8. — The distributions of bolometric corrections. The dotted histogram shows the corrections from JAi,3 Mm to L; r ; the dashed 
histogram is the corrections from vLs^ m to Lboli an d * ne solid histogram represents the bolometric corrections from v ^ 5100 ^ to ^bol as 
tabulated by R06. The legend lists the mean and standard deviation of the mean for each distribution. The overall bolometric corrections 
from fZ/3 Mm show a comparable fractional dispersion to the standard (5100A) bolometric corrections without being subject to dust extinction 
effects. 



